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ATLAS recently reported a Scr excess in a leptonic-Z + channel. This was interpreted in 

the literature in a simplified General Gauge Mediation model containing a gluino, a higgsino next- 
to-lightest supersymmetric particle (NLSP) and a gravitino lightest supersymmetric particle (LSP). 

We test the consistency of this explanation in lieu of the results of the corresponding search in 
GMS, and other LHG searches for New Physics. Due to non-decoupling effects from squarks the 
parameter space of these models is split into two regions; in one region additional leptons via top 
quark production is expected, while the other region sees a large probability for zero-lepton events. 

After combining the relevant constraints we find that these models cannot explain the ATLAS 
excess. 

PACS numbers: 12.60.Jy, 13.15.tg, 14.80.Ly 


I. INTRODUCTION 

A recent ATLAS search for beyond the standard model 
physics in a channel with two leptons, consistent with 
the production of a Z-boson, large missing transverse 
momentum and at least two jets, reports a 3 cr 

excess [1] for 20.3 fb“^ of integrated luminosity at a cen¬ 
ter of mass energy of 8 TeV. The other general purpose 
LHC experiment, CMS, has reported on a similar search, 
also with the full Run-I data set [2], seeing no excess. 
However, the cuts used in the two searches are different, 
and the observed ATLAS excess may a priori be consis¬ 
tent with the CMS results. Also, because the cuts are 
different, for some particular interpretation in terms of a 
new physics model one expects the predicted signal rates 
in each analysis to depend upon the signal kinematics. 
Hence, the relative number of predicted signal events in 
ATLAS as compared to CMS will depend in general upon 
the assumed interpretation, as well as on its parameters. 

In this article we investigate the consistency of the 
ATLAS excess with the CMS analysis, and with other 
searches at the LHC, for a General Gauge Mediation 
(GGM) model with a gravitino LSP [3]. The simplified 
model used here is inspired by [4] , and contains only three 
free parameters: the gaugino mass M 3 , fixing the gluino 
mass, tau/d, the ratio of the two higgs field vacuum ex¬ 
pectation values, and /i, the superpotential parameter, 
giving the mass of the higgsino NLSP, as well as one ad¬ 
ditional neutralino and one chargino, both dominantly 
higgsino. The excess can then be interpreted as stem¬ 
ming from gluino pair production and the decay chain 
g qqXi qqZG, depicted in Fig. 1 , with a leptonic 
Z-decay. In doing so we follow the model chosen by AT¬ 
LAS to interpret the results of their analysis. 

As the leptonic branching ratio of the Z is small, this 
model may also come in to conflict with recent searches 
for supersymmetry (SUSY) via jets and missing energy 
channels. This was briefly commented on in [5], but no 
detailed analysis of the parameter space was performed, 
and the article goes on to interpret the ATLAS results 
in an alternative model with Zs coming from the decay 


q q 



FIG. 1. Hypothesized GGM decay mode contributing to the 
ATLAS excess. 

X 2 ^ The ATLAS excess was also interpreted in 

[ 6 ], using two benchmark points in a GGM model with 
properties very similar to the model used by ATLAS, and 
in [7], in the context of a composite Higgs/Randall Sun- 
drum model, with heavy Kaluza-Klein gluon resonances 
decaying to vector-like quarks. 

II. MODEL 

For the model used in the ATLAS interpretation 
the remaining minimal supersymmetric standard model 
(MSSM) parameters were set as follows: the gaugino soft 
masses Mi = M 2 = 1 TeV, the sfermion soft masses 
= 1.5 TeV, and a gravitino mass light enough for the 
NLSP decays to be prompt. With rriq = 1.5 TeV, squark- 
gluino production dominates over gluino pair production 
for gluino masses above ^ 1 TeV. Squark-squark produc¬ 
tion would dominate sparticle production for high enough 
gluino masses. Also, with Mi = M 2 = 1 TeV more com¬ 
plicated decay chains open up for gluinos in this mass 
range, and the NLSP will no longer be dominantly hig¬ 
gsino for values of g close to 1 TeV.^ 


^ Alternatively, MSSM scenarios with a wino (or bino) NLSP 
could be considered, however, they have BR{xi —>■ jG) > 
0.23(0.77) [8], which should be easy to exclude from 
searches. 
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In order to explore a wider range of fi and M 3 val¬ 
ues without introducing such complications in the phe¬ 
nomenology we adopt a simpler model with Mi = 
M 2 = 1.5 TeV and sfermions completely decoupled at 
mj = 4.5 TeV, keeping in mind that lowering the squark 
mass scale generally will lead to stronger bounds on the 
model. The mass parameters are defined at a scale of 
~ 4.5 TeV. 

The gravitino mass is given by the scale of SUSY 
breaking, but must be very light for the NLSP to de¬ 
cay promptly; we set it to be effectively zero for the col¬ 
lider simulation. All results are presented for tan/3 = 1.5 
and tan/3 = 30. The choice of low tan/3 is made in or¬ 
der to maximize the branching ratio Xi ZG^ which 
is in competition with ^ f® a smaller ex¬ 

tent ^ 7 G. For low values of tan/3, /i > 0 and a 
higgsino NLSP, this is approximately 100% [8-10]. In¬ 
creasing tan^d will decrease the signal. The two values 
used thus explore different parts of the parameter space. 

The lightest higgs mass is simply set to the experi¬ 
mentally measured value of rrih = 125.09 ±0.24 GeV [ 11 ] 
by assuming extra operators in the higgs sector, e.g. by 
using dimension-5 operators as proposed in [12]. 

The relative squark masses (and to a smaller extent 
the value of Mi and M 2 ) determine the branching ra¬ 
tio of the gluino decay into the various quark flavours. 
Since the NLSP is dominantly higgsino there will nec¬ 
essarily be large branching ratios into third generation 
quarks. When these are kinematically forbidden, the loop 
induced decays g gxi 2 become important. Due to the 
importance of decays involving third generation quarks, 
tan/3 also affects the gluino branching ratios. Further 
complicating matters is the existence of multiple higgsi- 
nos at roughly the same mass [xi^ X 2 XiY • In Fig. 2 
we show the branching ratios of the gluino calculated us¬ 
ing SUSYHIT 1.4 [15], as a function of Am = rrig — m^o 
for tan/3 = 1.5 (solid lines) and tan= 30 (dashed lines). 
The gluino mass is fixed at rUg = 900 GeV and the other 
parameters are as given above. 

We see that at least one top quark will be produced per 
event on average even down to Am ^ 350 GeV, where 
the proximity of the g thxf threshold becomes impor¬ 
tant. For tau/d = 30 we get a sizeable contribution from 
the decays g 56x5 2 5 reducing somewhat the produc¬ 
tion of top quarks at high Am. We find that iowering 
Ml and/or M 2 down to the giuino mass changes iittie: 
there is a siight increase in the first and second generation 
quark decays versus the giuino ioop decay beiow the thxf 
threshoid, but no significant impact above Am ^ 350 
GeV. This ensures the presence of a significant number 
of events with additionai ieptons from ieptonic top de¬ 
cays, which we wiii see have an impact on the aiiowed 


^ We note that a co-NLSP scenario is unlikely for higgsinos. This 
would require relatively small but negative Mi values and large 
M 2 , see [13, 14]. 



Am [GeV] 

FIG. 2. Branching ratios for the gluino as a function of 
the gluino-neutralino mass difference with tan/3 = 1.5 (solid 
lines) tan/3 = 30 (dashed lines). The gluino mass is fixed at 
Trig — 900 GeV. Some lines are not visible because they are 
at very low values of the gluino branching ratio. 

parameter space of the model. The decays to first and 
second generation quarks are heavily suppressed. 

The ATLAS analysis assumed equal branching frac¬ 
tions of ^ ^ for Q = ignoring the heavy 

quark decays. This simplifying assumption has relatively 
little impact on their analysis and the bounds set be¬ 
cause of the focus on leptons from the Z boson. However, 
the structure of GGM predicts generic sum rules on the 
sfermion soft masses [3], 

rriQ — 2m}j ± m% — m\ ± m% = 0 (1) 

2rriQ — mfj — mjj — 2m\ ± m\ = 0, (2) 

which makes decoupling only the third generation 
squarks challenging.^ One is faced with the choice of ei¬ 
ther making the phenomenological assumption that the 
sum rules are broken somehow, increasing the number of 
free soft sfermion mass parameters in the model, and 
requiring an explanation of why the third generation 
masses are significantly heavier, e.g. one could speculate 
that something along the lines of higgsed gauge media¬ 
tion [16] could work.^ Or, to reduce the number of top 
quarks produced, one can at best decouple all squarks ex¬ 
cept the lightest bottom squarks by an appropriate choice 


^ Technically the ATLAS model does not fulfil these sum rules, 
however, a slight modification of the soft mass parameters would. 
With first and second generation sfermions at 1.5 TeV, the third 
generation sfermions must be raised to ~ 5 TeV for the gluino 
branching ratios to light quarks to equal those to third generation 
quarks. 
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observed 29 

background 10.6zb3.2 

number of sigma 3.0 

s (95% CL) 7.1-31.8 

TABLE I. Summary of ATLAS_onZ constraints, showing the 
observed number of events, the number of expected Standard 
Model events inferred from data, the number of sigma the 
excess corresponds to and the 95% CL constraint upon a pu¬ 
tative number of signal events s. The first three data are taken 
from Ref. [1], whereas we infer the bound on s ourselves (see 
text). 

of the soft masses, however, the decay g tbxf remains. 
Neither seems very consistent with the simple model of 
ATLAS. As a result we will here include gluino decays 
to top and bottom quarks in the model. 


III. SCAN AND SIMULATIONS 

We perform a grid scan over the range 0-1500 GeV in 
g and M 3 , using a step size of 15 GeV in both directions. 
At each step we calculate the resulting sparticle spectrum 
using SOFTSUSY 3.5.1 [17] and the sparticle branching 
ratios with SUSYHIT 1.4 [15]. Spectrum and decay in¬ 
formation is communicated via the SUSY Les Houches 
Accord [18], using PySLHA [19]. For all parameter points 
we check that rrig > m^o, and that the NLSP is mostly 
higgsino (more than 0.90). At each point we generate 
100 000 SUSY Monte Garlo events with gluino pair pro¬ 
duction using Pythia 8.186 [20, 21]. The cross sections 
used are based on Prospino [22], using the NLLfast soft¬ 
ware including also NLL re-summation of soft gluon emis¬ 
sion [23-26]. These events are then propagated through 
our implementations of several collider analyses, detailed 
below. 

The ATLAS analysis with the excess requires two lead¬ 
ing opposite sign same flavor (OSSF) leptons with pt > 
25,10 GeV and invariant mass 81 < mu < 101 GeV, a 
minimum missing transverse energy of > 225 GeV, 
at least two jets, and total transverse energy Ht > 
600 GeV, where Ht is given as the scalar sum of the 
transverse momenta of the two leading leptons and all 
accepted jets. Jets are reconstructed with the anti-/cT 
algorithm [27] using Fast Jet [28], with a jet radius pa¬ 
rameter of R = 0.4, and are required to have pr > 35 
GeV and lie within I 77 I < 2.5. In the following we will de¬ 
note the signal region with the sum of ee and pp events 
as ATLAS_onZ. We summarize the ATLAS measurements 
in Table 1. In order to calculate a constraint upon the 
number of non-Standard Model signal events s, we pro¬ 
file over a Gaussian background rate, see Section IV, but 
otherwise use Poisson statistics. If other combined con¬ 
straints predict an ATLAS_onZ signal rate outside of this 
range, we shall conclude that they are incompatible with 
the signal at the 95% GL. 


The on-Z GMS analysis [2], here called CMS_onZ, re¬ 
quires a leading pair of opposite sign same flavor lep¬ 
tons satisfying pt > 20 GeV and 81 < mu < 101 GeV. 
Three signal regions are constructed, covering the ranges 
100-200 GeV, 200-300 GeV and > 300 GeV in 
all requiring at least three jets with pt > 40 GeV and 
\rj\ < 3.0. For jet reconstruction the anti-/cT algorithm 
with R = 0.5 is used. A notable difference with respect 
to the event selection in ATLAS_onZ is that no cut on Ht 
is applied. While there are some details of the original 
analysis which are difficult to reproduce outside of the 
experimental collaborations, e.g. trigger efficiencies, by 
simulating models similar to those used for interpreta¬ 
tion in ATLAS_onZ and CMS_onZ we have checked that our 
implementations reproduce the observed limits to within 
theoretical uncertainties, under the assumptions made. 

In addition to the leptonic-F + F^^®® analyses from 
GMS and ATLAS, the scenario used here could be con¬ 
strained by other searches involving leptons. This in¬ 
cludes three and four lepton final states where extra lep¬ 
tons are produced in leptonic top decays, or from two 
chains with Fs. The latter is heavily suppressed by the 
leptonic branching ratio of the F, down to ^ 7% of the 
number of events with a single leptonically decaying F, 
thus of the order of two events could be expected for the 
given luminosity, depending on the exact cuts of such an 
analysis. 

We check the most relevant searches which are the 
ATLAS stop search with leptons ATLAS_stop_L100 [29], 
and the GMS multi-lepton search with three or four 
leptons, CMSjnultilepton [30]. From the analysis in 
ATLAS_stop_L100 we include the signal region LlOO re¬ 
quiring exactly two opposite-sign leptons with pt > 
25,10 GeV, at least two jets with pt > 100, 50 GeV, a 
‘stransverse mass’ mT 2 >100 GeV and an invariant mass 
mu for the two leptons outside the range 71-111 GeV. 
For the GGM model studied here, the cut on mu means 
that ATLAS_stop_L 100 is mainly sensitive to events where 
neither of the two Fs decay leptonically. 

The GMS multi-lepton search requires at least three 
isolated leptons with pt > 20,10,10 GeV within |? 7 | < 
2.4. Jets are subject to the requirements pt > 30 GeV 
and 1 7 ^I < 2.5. Accepted events are divided into a 
large number of signal regions based on the number of 
opposite-sign same-flavour lepton pairs, F^^®®, the pres¬ 
ence of a OSSF lepton pair with an invariant mass in the 
75-105 GeV range, the scalar sum of jet pts and the num¬ 
ber of tagged 6-jets. Due to an overlap between CMS_onZ 
and the most relevant signal regions in CMS jnulti lepton, 
in the combination we use CMS jnulti lepton for param¬ 
eter regions where the gluino-neutralino mass difference 
Am is larger than 500 GeV and CMS_onZ for mass differ¬ 
ences smaller than this. A choice like this is necessary 
in order to have statistically independent signal regions 
when we do not have the information to take into account 
the correlations. As can be seen from Fig. 2, the choice 
ensures that CMS jnulti lepton is only applied in the re¬ 
gion of parameter space where additional leptons can be 
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expected due to a sizeable production of top quarks. 

Given the small leptonic branching ratios of W and Z, 
the GGM scenario studied can also be constrained from 
searches for zero-lepton final states. We include the sig¬ 
nal region 3j from the ATLAS search for final states with 
jets and missing energy [31], here called ATLAS_jMET_3j. 
Besides a lepton veto, this signal region requires at least 
three jets with pt > 130, 60, 60 GeV and \r]\ < 2 . 8 , miss¬ 
ing energy > 160 GeV and an ‘effective mass’ 

rueff > 2200 GeV. Also, the missing energy is required 
to account for at least 30% of the effective mass combi¬ 
nation of and the three leading jet pts. 

No significant excesses are seen in either of these 
searches. For all three searches we have checked, as 
above, that we can reproduce the relevant limits on SUSY 
interpretations presented by the experiments. 


IV. STATISTICS 

In order to combine the results from all the analy¬ 
ses, each independent signal region i is assigned a likeli¬ 
hood Ci consisting of a Poisson factor for the total event 
count and a Gaussian for modelling the background un¬ 
certainty: 


= Pois(nj|sj + bi) x Gauss(&„J&j,(TbJ. (3) 

Here rii is the observed number of events, Si and hi are the 
expected number of signal and background events, and 
hrrii is the observed background measurement with an 
expected standard deviation . Inserting the observed 
values for n^, 6 ^. and cr^. we are left with a likelihood 
function for the two parameters Si and hi. While Si will 
be a function of the SUSY parameters p and M 3 , hi is 
an unknown nuisance parameter which we eliminate by 
profiling Ci over hi. With all the Ci coming from inde¬ 
pendent signal regions, the combined likelihood is then 
simply given by 


£(s)=(4) 

i 


where the double hat indicates that we have maximized 
^i{si^bi) over hi subject to a fixed value of Si. 

For any given parameter point in p and M 3 , the signal 
expectation values s are in principle fully determined. 
In order to set limits in the model parameter space we 
introduce a common signal strength parameter ps such 
that the expected signal yield in signal region i is PsSi. 
Points in the SUSY parameter space for which the upper 
limit on ps is found to be less than 1 will be excluded at 
the confidence level chosen for the test. 

For every choice of the SUSY parameters we now have 
a single-parameter likelihood function C{iis) = >C(/isS). 
From the likelihood ratio 




^(As)’ 


(5) 


we construct a test statistic q given by 


-2\nX{ps) 

0 


As ^ Ms 
As ^ Ms 5 


( 6 ) 


where As is the value of ps that maximizes C{fis), i-e. 
the signal strength value preferred by the observed data. 
Higher values of q correspond to increasing disagreement 
between data and the hypothesized value of /i^, but only 
in the direction of /ig > As- For a given ps the observed 
value g'obs of Q is calculated from the data. The p-value 
for this observation is then found from 

00 

= J f{Q\P's)dq, (7) 

Qohs 


where f{q\jiis) is the pdf of q. To determine p^^ we make 
use of the asymptotic limit in which /(gj/ig) is given by 
a ‘half chi-square’ distribution, i.e. an equally-weighted 
sum of a delta function at zero and a chi-square distri¬ 
bution for one degree of freedom [32] . The 95% GL^ [33] 
upper limit on ps is the highest value of ps satisfying 

> 0.05, (8) 

1 -Po 

where po is the p-value for the test statistic 


j-21nA(0) /is>0 

^0 ^ 1 n ' ^ n 

1^0 ^ 0, 

used to test the level of disagreement between data and 
the background-only hypothesis. 


V. RESULTS 

We show the 95% GL allowed region for the 3cr 
ATLAS_onZ excess in Fig. 3 for tan A = 1.5,30 as the 
lighter band. The 95% CL^ excluded regions from the 
other searches are overlaid, and for reference the two AT¬ 
LAS benchmark points at = (700,200) GeV 

and (900,600) GeV, tan A = 1-5, are indicated with 
white diamond markers. It is clear from the figure 
that all points explaining the ATLAS_onZ excess for ei¬ 
ther value of tan A fall afoul of at least one of the other 
searches. Indeed, CMS_onZ alone is already incompati¬ 
ble with the excess at the 95% GL^ level except for two 
small regions. One region, around m^o « 950 GeV, 
nig ~ 980 GeV which is anyway well excluded by the 
ATLAS-j MET searches. The other region, including the 
point m^o = 190 GeV, nig = 930 GeV is excluded by 
CMS jnultilepton. When tan A is changed, the exclu¬ 
sion contours move somewhat in gluino and lightest neu- 
tralino mass, but the qualitative conclusions remain un¬ 
changed. 

Given the tension between the other searches and the 
ATLAS_onZ excess, we now combine the CMS_onZ exclu¬ 
sion with various other searches to see what the combined 
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ATLAS scenario. Still, both ATLAS benchmark points 
are excluded at the 95% confidence level from the com¬ 
bination of ATLAS_onZ and CMS_onZ alone. On the other 
hand, there are still some points left allowed at the 95% 
CL that predict an ATLAS_onZ signal rate of up to 13(12) 
for tan/3 = 1.5(30). These are within the 95% CL signal 
rate region of 7.1 — 31.8 and so are still compatible with 
the ATLAS_onZ signal at the 95% CL. 



600 -^^^^- Mq 

200 400 600 800 1000 1200 

m^o [GeV] 


FIG. 4. The 95% CL^ exclusion curves in the {m^o,mg)- 
plane for tan/3 = 1.5 (white) and tan/3 = 30 (black), using 
both the ATLAS_onZ and CMS_onZ signal regions. The color 
map shows the expected number of ATLAS_onZ signal events 
for the tan/3 =1.5 scenario. The two ATLAS benchmark 
points are indicated with white diamond markers. The region 
below each curve is excluded by the combination. ATLAS_onZ 
signal events are constrained to be below 13(12) at the 95% 
CL for tan/3 = 1.5(30). 


FIG. 3. The band in the (m^o,m^)-plane preferred 
by ATLAS_onZ at 95% CL (lighter region), compared to 
the coloured 95% CL^ exclusion contours from CMS_onZ, 
ATLAS_stop_L100, CMSnnultilepton and ATLAS_jMET for 
tan yd = 1.5 (top), tan/3 = 30 (bottom). The ATLAS_stop_L100 
exclusion region boundary is shown as a dashed line. Two 
ATLAS benchmark points are indicated with white diamond 
markers. 


data set predicts for the number of signal events in the ex¬ 
cess. The 95% CLg bound in the (m^o, m^)-plane result¬ 
ing from combining ATLAS_onZ and CMS_onZ is shown in 
Fig. 4. Here the white contour depicts the limit obtained 
for tan P = 1.5, while the black contour is for tan /3 = 30. 
Also shown are the bounds given a 20% systematic un¬ 
certainty on the cross section. The color map shows the 
predicted number of signal events for the ATLAS_onZ anal¬ 
ysis in the scenario with tan/3 = 1.5. Since the squarks 
are decoupled from gluino production here, the produc¬ 
tion cross section for the model is lower than the original 


If we add on the contributions from ATLAS_stop_L100 
and CMSjnultilepton the resulting 95% CLg bound is 
shown in Fig. 5. As expected, the exclusion limits are 
improved in the regions of large gluino-neutralino mass 
difference, where the production of additional leptons 
through top quarks is significant. The slight dip in the 
contour at (m^,m^o) ^ (1000,500) GeV is where the 
domains of the CMS jnulti lepton and CMS_onZ analy¬ 
ses meet. In the region of small gluino-neutralino mass 
difference the limit remains approximately unchanged, 
and the combined allowed region predicts up to 13(11) 
ATLAS_onZ signal events for tan/3 = 1.5(30), still consis¬ 
tent with the 7.1 — 31.8 95% CL constraint. 

The final exclusion limit, obtained after including also 
ATLAS_jMET_3j , is shown in Fig. 6. Due to the lepton veto 
in this analysis, the exclusion limit is mainly strength¬ 
ened in the region with gluino-neutralino mass differ¬ 
ences less than 400 GeV, where the main source of lep¬ 
tons is through the small leptonic branching ratio of the 
Zs. With tan yd = 30 this effect is further enhanced by 
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the reduced branching ratio into Zs. We note that gluino 
masses below 1 TeV are fully excluded for both values of 
tan/3, and the remaining allowed parameter space has 
a maximum of 6 expected signal events for ATLAS_onZ, 
which is far from explaining the observed excess. 



FIG. 5. The 95% CLs exclusion curves in the (m^o,m^)- 
plane from combining ATLAS_onZ, CMS_onZ, ATLAS_stop_L100 
and CMSjnultilepton. The color map shows the expected 
number of ATLAS_onZ signal events for the tan/3 = 1.5 sce¬ 
nario. The two ATLAS benchmark points are indicated with 
white diamond markers. The region below each curve is ex¬ 
cluded by the combination. ATLAS_onZ signal events are con¬ 
strained to be below 13(11) at the 95% CL for tan j3 = 1.5(30). 

The observed excess in ATLAS_onZ was also recently in¬ 
terpreted within a GGM framework in [6], where two new 
benchmark points are presented. The main difference 
with respect to the ATLAS scenario is heavier squarks, 
and that the lightest neutralino is a wino-bino mixture. 
The first point, referred to as GGMl, has rrig = 1088 GeV, 
m^o = 428 GeV and squark masses around 2800 GeV. 
The second point, GGM2, has a higher production cross 
section due to lighter gluinos and squarks, at 911 GeV 
and ^ 2400 GeV, respectively. Combining the collider 
constraints considered here we find that GGM2 is excluded 
at the 95% confidence level, while GGMl escapes exclusion. 
However, it should be noted that GGMl predicts only ^ 3 
signal events for ATLAS_onZ^. 

One may ask: can one tweak the simplified model in 
order to squeeze around the constraints? In our anal¬ 
ysis, we have set the simplified model up in order to 
maximize the ATLAS_onZ region compared to the other 
constraining searches. If tan/3 is increased further, the 


^ We note that this is in slight disagreement with the simulation 
in [6], which finds 6±1 predicted signal events for GGMl. However, 
this discrepancy does not change the conclusion. 



m^o [GeV] 

FIG. 6. The 95% GL^ exclusion curves in the (m~o,m^)- 
plane from combining all collider searches detailed in the text. 
The color map shows the expected number of ATLAS_onZ sig¬ 
nal events for the tan/3 = 1.5 scenario. The two ATLAS 
benchmark points are indicated with white diamond markers. 
The region below each curve is excluded by the combination. 
ATLAS_onZ signal events are constrained to be below 6(5) at 
the 95% GL for tan/3 = 1.5(30). 


ATLAS_onZ signal decreases for the same gluino/lightest 
neutralino masses. Thus, these would have to be lowered 
in order to get a signal to fit, and such lighter sparti- 
cles would suffer more from the other searches. If we 
were to make squarks lighter, although the ATLAS_onZ 
signal would increase for the same gluino/lightest neu- 
tralino masses, the ATLAS_jMET constraints would become 
much stronger. In any case, CMS_onZ is in tension with 
nearly all of the ATLAS_onZ parameter space, and this is 
unlikely to change. The CMSjnultilepton constraint is 
mainly due to events with one leptonically decaying Z- 
boson plus additional leptons from decaying top quarks, 
predicted by the signal model due to the higgsino nature 
of the neutralino. If the model is manipulated to reduce 
the production of such additional leptons, the weaken¬ 
ing of the CMS jnulti lepton constraint will be compen¬ 
sated by a corresponding strengthening of ATLAS_j MET as 
more events will pass the lepton veto. Thus, although we 
have not exhaustively covered the full multi-dimensional 
MSSM plus light gravitino space, there are good grounds 
for expecting that our conclusions — that the ATLAS_onZ 
excess is incompatible with other searches at the 95% CL 
for GGM-type models — also apply to the full MSSM 
space. 
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VI. CONCLUSIONS 


In conclusion, we have seen that a simplified GGM 
model with only M 3 and /i as free parameters, and with 
higgsino y?, X 2 ^ gravitino LSP as the 

only sparticles produced at the LHG, cannot explain the 
ATLAS excess reported in [ 1 ] when faced with results 
from other current new physics searches. Strong bounds 
on the model that can be set from other leptonic searches 
are due to the higgsino nature of the NLSP, leading to 
the production of top quarks with leptonic decays. Zero- 
lepton searches also provide strong constraints in the 
parameter regions where leptons are mainly produced 
through the small leptonic branching ratio of the Z. Ten¬ 
sion between ATLAS _onZ and the other relevant searches 
is evident in Fig. 3. A combined fit to all constraints. 


including the excess^ at the 95% GL predicts less than 6 
signal events in the ATLAS_onZ search region, compared 
to 7.1-31.8 being inferred from the ATLAS_onZ search re¬ 
gion alone. 
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